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ABSTRACT 

We discuss the origin of galactic angular momentum, and the statistics of the present 
day spin distribution. It is expected that the galaxy spin axes are correlated with the 
intermediate principal axis of the gravitational shear tensor. This allows one to re- 
construct the shear field and thereby the full gravitational potential from the observed 
galaxy spin fields. We use the direction of the angular momentum vector without any 
information of its magnitude, which requires a measurement of the position angle and 
inclination on the sky of each disk galaxy. We present the maximum likelihood shear 
inversion procedure, which involves a constrained linear minimization. The theory is 
tested against numerical simulations. We find the correlation strength of nonlinear 
structures with the initial shear field, and show that accurate large scale density recon- 
structions are possible at the expected noise level. 

Subject headings: galaxies: statistics — large-scale structure of universe 



1. Introduction 

In the gravitational instability picture of structure formation, galaxy angular momentum arises 
from the tidal torquing during the early protogalactic stages (Hoyle 1949). Gravitational torquing 
results whenever the inertia tensor of an object is misaligned with the local gravitational shear 
tensor. One might thus expect that the observed galaxy spin field would contain some information 
about the gravitational shear field, and possibly allow a statistical reconstruction thereof. Although 
the observational search has yet to detect any significant alignments of the galaxy spins so far 
(Han et al. 1995), it has been shown that the sample size has to be increased in order to detect 
weak alignments (Cabanela & Dickey 1999 and references therein). 
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It was argued that angular momentum only arises at second order in perturbation theory 
(Peebles 1969). This picture involved a spherical protogalaxy, which at first order has no quadrupole 
moment, and thus cannot be tidally torqued. Later, it was realized that protogalaxies should be ex- 
pected to have order unity variations from spherical symmetry, and should be tidally torqued at first 
order (Doroshkevich 1970; White 1984). However, recent quantitative predictions of this picture 
(Catelan & Theuns 1996, hereafter CT) appear to be contradicted by simulations (Lemson & Kauffmann 1999), 
which leaves the field in a state of confusion. 

The Eulerian angular momentum of a halo in a region Ve relative to its center of mass in 
comoving coordinates, Jj Eul = J y ^ px x v<i 3 x, can be written in terms of Lagrangian variables to 
second order accuracy using the Zel'dovich approximation: 

L 2nd oc - f pq X V^d 3 q, (1) 

Jv L 

where Vl is a Lagrangian counter part of Ve, and 4> is the gravitational potential. Approximating 
in equation (1) by the first three terms of the Taylor-series expansion about the center of mass, 
we obtain the first order expression of the angular momentum (White 1984) in terms of the shear 
tensor T = (Tji) = (djdi<p) and the inertia tensor I = (_%) = (f qiq^q): 

Lj st oc e ijk TjiIi k . (2) 

CT used equation (2) to calculate the linearly predicted angular momentum under the approxima- 
tion that the principal axes of I and T are uncorrelated. They also discussed the small factor by 
which the neglect of the correlation between I and T overestimates the angular momentum in the 
context of the Gaussian-peak formalism. In § 2, we show by numerical simulations that this factor 
in fact is dominant. 

For galaxies, the direction of angular momentum is relatively straightforward to observe, while 
the magnitude very difficult. We will thus concentrate on the statistics of the direction of the spin, 
dealing with unit spin vectors. We note that the spin in equation (2) does not depend on the trace 
of either I nor T. So we can consider unit trace-free tensors, (TijTij = 1), and the one point 
statistics of the spins. The most general quadratic relation between a unit spin vector and a unit 
traceless shear tensor is ^ 

(LjLjjT) = Qij = — - — Sij — aTikTkj, (3) 

where a 6 [0, 1] is the correlation parameter, measuring how well aligned I and T are. If I and T 
are uncorrelated, a = 1. If they are perfectly correlated, a = 0. It is also a convenient parameter- 
ization of the nonlinear effects on the spin-shear correlation, and must be measured in numerical 
simulations. 
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2. Simulations 

We ran the N-body simulations using the PM (Particle-Mesh) code described by Klypin & 
Holtzman (1997) with 128 3 particles on a 256 3 mesh in a periodic box of size Lb = 80/j _1 Mpc 
for the cold dark matter model (Qq = 1, ag = 0.55 and h = 0.5). We identified halos in the 
final output of N-body runs (total number of halos, N t = 2975) by the standard friends-of-friends 
algorithm with a linking length of 0.2. The angular momentum of each halo was measured from the 
positions and velocities of the component particles in the center-of-mass frame. We also calculated 
the initial shear tensor by taking a second derivative of the gravitational potential at the Lagrangian 
center-of-mass of each halo. 

Through the simulations, we first tested the validity of the linear perturbation theory by 
calculating the average correlations among the 1st, the 2nd order Lagrangian and the Eulerian 
angular momentum obtained from simulations. We have found that (L 2nd • L Eul ) = 0.55, (L lst • 
L Eul ) = 0.51, and (L lst • L 2nd ) = 0.62. It shows that nonlinear effects only add a factor of 2 scatter 
to the linearly predicted spin-shear correlation (Sugerman et al. 1999). 

Rotating the frame into the shear-principal axes, we find the optimal estimation formula for 
the parameter a from equation (3): 

3 

a = 2-6^A?|L;| 2 , (4) 
i=i 

where {Aj} 3 =1 are the three eigenvalues of the trace-free unit shear tensor, satisfying A« = 1 
and \ = 0- Using equation (4) along with the initial shear tensors and the angular momentum 
of final dark halos measured from the simulations, we calculated the average value of a to be 
0.237 ± 0.023, showing 10<r a deviation [a a = y/8/ (5-/V t )] from the value of 0. The inferred value 
of a results in a significant correlation of the shear intermediate principal axis with the direction 
of the angular momentum of halos. We have shown by simulations that the shear and the inertia 
tensors are misaligned with each other to a detectable degree, generating net tidal torques even 
at first order, although they are quite strongly correlated (it was found by the simulations that 
(In) = —0.85, (I22) = 0.11, (133) = 0.74 for the trace free unit inertia tensor I in the shear principal 
axes frame). Thus we conclude that the CT approximations are useful for determining the direction 
of the spin vector, but poor for predicting its magnitude (Lemson & Kauffmann 1999) due to the 
strong correlation between I and T. The strong correlation can be understood physically if one 
considers I to be the collection of particles which have shell crossed, in which case it is identical to 
T. 
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3. Shear Reconstruction 

As in the CT prescription, the probability distribution of P(L|T) can be described as Gaussian, 
and its directional part, P(L|T) is calculated by 

P(L|T) = J P(L\T)L 2 dL = IQ |J /2 (i'Q^ty 3 ^ (5) 

where Q is given from equation (3). According to Bayes' theorem, P(T|L) = P(L|T)P(T)/P(L). 
Here, P(T) = P[T(xi), T(x2), . . .] is a joint random processes linking different points with each 
other. In the standard picture of galaxy formation from random Gaussian fields, P(T) is given as 
P(T) = Aexp(-T*C~ 1 T/2) where C = (C ijk i) = (T ij (x)r H (x + r)) is the two-point covariance 
matrix of shear tensors: 



(6) 



Cijki{r) = (SijSki + SikSji + 5ii5 jk )^Y ~ y^j + {^/kn)^(r) + ~Y -~Y 



Here f = r/r, J n = nr n £(r)r n 1 dr, and ^(r) = (<5(x)<5(x + r)) is the density correlation 
function. 

We would like to compute the posterior expectation value, for example (T|L). By symmetry, 
both the expectation value and maximum likelihood of T occur at T = 0, which is not the solution 
we are looking for. We must thus consider the constrained expectation, with the constraint that 
J TijTijd 3 *. = 1. In the limit that P(T|L) is Gaussian, this is given by the maximum eigenvector 
of the posterior correlation function (jj; m (x a ,x^) = (Tjj(x a )r; m (x i g)|L) . The solution satisfies 



/ 



^m(x a ,X ( g)T ii (x a )(i 3 X a = AT im (x /3 ) (7) 



where A is the largest eigenvalue for which equation (7) holds. In the asymptotic case that a ^ 1 as 
in our simulations results, we have a simple expression for the LjLj-dep ending part of the posterior 
correlation function in terms of the traceless shear, Ty = Ty — 5ijTu/3: 

iijlm( x a, x f3) = — J Cij n k(X- a — X^)Cl mo k(^-/3 ~ X 7 )L n (x 7 )L (x 7 )(i 3 X 7 , (8) 

where C is now the trace- free two-point covariance matrix of the shear tensors, related to C by 
Cijki = Cijki-SkiCi jnn /3-5ijC mm ki/3 + 5ij5kiC mmnn /9. Substituting (3) into (8) explicitly satisfies 
a modified version of (7). It is readily described in Fourier space: 



/ 



lM^Mf lj{K)Aa = AXUV). (9) 



Equation (9) differs from (7) by the appropriate Wiener filter due to small a, since we have assumed 
no noise. 
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Since we have discarded the constant diagonal component which does not depend on LiLj, 
equation (8) is not positive definite, and in fact has zero trace. We can nevertheless use a power iter- 
ation to quickly obtain the eigenvector corresponding to the largest eigenvalue. One starts with an 

initial guess Tg(x Q ), and defines iterates f™ +1/2 = J ^, m 2^, and f^ 1 = f^ +1/2 /^J(f^ +1/2 y + 
T^" 1 . This effectively eliminates the negative eigenvectors from the iteration, and converge to the 
correct solution. After n iterations, the fractional error is proportional to (Ai/Ao) n , where Ao and 
Ai are the largest and second largest eigenvalues, respectively. In our experience, fifty iterations 
converges the result to within a percent of the exact solution. In practice, one can expect the corre- 
lation to drop off rapidly at large separations, requiring only an evaluation of a sparse matrix- vector 
multiplication. 

The final step is the reconstruction of the density given the traceless shears. It is convenient 
to consider the full shear as a orthonormal reparametrization in terms of trace and trace-free 
components into a vector Tj = {S/y/3, V2, V3, V2T23, \f2T%i\ where 5 = T\\ + T22 + 

v 2 = [(-3 - V3)T U + 2V3T22 + (3 - V3)T 33 ]/6 and v 3 = [(3 - y/3)T u + 2^3T 22 + (-3 - V3)T 33 ]/6. 
The two point correlation function (6) is just linearly transformed into these new variables, which 
we will denote (Tj(x)Tj(x + r)) = C'-(r), where the indices i,j run from 1 to 6. The inverse 
correlation will be denoted D = (C')" 1 . From D we extract a scalar correlation Du(r) and a 5 
component vector correlation = {-D12, -D13, D14, D15, D\q}. Similarly, we will use a Greek index 
to denote the 5 component vector of Tj where i > 1. We then obtain the expression for the 
density (5\T^) = D^E^T^. If the galaxies are uniformly spaced on a lattice, one could also use a 
Fast Fourier Transform to rapidly perform the shear reconstruction iterations, and the projection 
of traceless shear to density. 

To test the reconstruction procedure, we have made a large stochastic spin field on a 128 3 
lattice and applied the reconstruction procedure with a Wiener filtering scale of k = 32 for 
the case of a power-law spectrum, P(k) = k~ 2 . Figure 1 plots the correlation coefficient c = 
(<5(k)<5 T -( — k))/-y/ P(k)(\5 r (k)\ 2 ) between the reconstructed density field S r and the original density 
field 5. We see that for a = 1, corresponding to the case of uncorrelated shear and inertia tensor, 
the reconstruction is very accurate on large scale (small k) and noisy at small scale (large k) as 
expected. Even for the realistic case of a ~ 0.2 given the nonlinear effects, Figure 1 shows that 
a reconstruction is still quite possible but with less accuracy. Note that the reconstructed signal 
automatically becomes Wiener filtered on large scale. 

4. Observational Effects 

Let us now mention several complications that might arise with real data, and speculations on 
how they might be addressed. 

The three dimensional spin axis of a disk galaxy can in principle be determined (Han et al. 1995), 
but it requires knowledge of the rotation curve as well as morphological information, such as ex- 



-6- 



density reconstruction 
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Fig. 1. — Cross correlation coefficient c = (5(k)5 r (— k))/y / P(/c)(|(5 r (k)| 2 ) between original S and 
reconstructed S r density field in a 128 3 spin field simulation. The correlation parameter a = 1 refers 
to the purely random shear model (Catelan & Theuns 1996), while a = 0.2 to the realistic shear 
model as found in our N-body simulations. The large scale density field is always well reconstructed, 
while noise dominates on small scale. 
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tinction maps or the direction of spiral arms. A much easier observational task is to measure only 
the position angle a on the sky and the projected axis ratio r. Assuming a flat disk geometry 
for the light and the z coordinate along the radial direction, this corresponds to observation of 
r = \L Z \ and tan(a) = L y /L x . We immediately note two discrete degeneracies: L z — > —L z 
and {L x ,L y } — > —{L x ,L y }, giving four possible spin orientations consistent with the observ- 
ables. Since we only care about the spatial orientation of the spin axis, but not its sign, only a 
twofold degeneracy exists, which we call L a , L 6 . We can readily take that into account. Noting 
that P(L a V L b \T) = P(L a |T) + P(L 6 |T) - P(L a A L 6 |T), the last term is zero, we simply replace 
all occurrences of Lj-L,- in equation (8) with L^L - + L b L b . 

The reconstruction itself is quite noisy for an individual galaxy, so we will treat as noise 
the displacement between the observable galaxy positions in Eulerian redshift space and their 
Lagrangian positions used in our analysis. If galaxies are randomly displaced from their origin by 
ay, we simply convolve the two point density correlation function £(r) with a Gaussian of variance 
of, and appropriately update the shear correlation matrix. We furthermore note that the use of 
redshift for distance introduces an anisotropy in the smoothing Gaussian, which can also be readily 
taken into account. This procedure could in principle be inverted to measure the nonlinear length 
scale by measuring ay if one knows the intrinsic power spectrum. One could vary ay until one 
maximized the eigenvalue A in Equation (7). 

The shear reconstruction correlator (8) is only accurately defined at each galaxy position, so 
the integral is should be replaced by a sum over galaxies. The eigenvector iteration scheme works 
the same way with the discretized matrix. The density is recovered from the trace-free shear in the 
same fashion as the continuum limit, where again we only reconstruct the density at the galaxy 
positions. The density at any other point is reconstructed in a similar fashion using the posterior 
expectation value given the Gaussian random field. 

The observed spin vectors are of the luminous materials. They can be observed out to large 
radii, tens of kpc, through radio emission of the gas. While some warping of the disks are seen, the 
direction of the spin vector generally changes only very modestly as one moves to larger radii. This 
suggests that the luminous angular momentum be well correlated with the angular momentum of 
the whole halo. The standard galaxy biasing is not expected to affect the reconstruction, since we 
have not used the galaxy densities but only the spin directions. Merging of galaxies is also taken 
into account by this reconstruction procedure, since the resulting spin vectors are expected to align 
with the constituent orbital angular momentum vectors. The latter is again predicted by the same 
shear formula. The merger effects are included in the N-body simulations presented above. 

5. Conclusions 

We have presented a direct inversion procedure to reconstruct the gravitational shear and 
density fields from the observable unit galaxy spin field. The procedure is algebraically unique, 



-8- 



dealing only with linear algebra, and is computationally tractable when done iteratively. We have 
shown that the angular momentum-shear correlation signal in simulations is strong enough to allow 
useful reconstructions. Direct N-body simulations suggest that nonlinear collapse effects only reduce 
the linearly predicted spin-shear correlation by a factor of 2. The deprojection degeneracy can be 
incorporated in a straightforward fashion. Large galaxy surveys with a million galaxy redshifts 
are coming on line soon. This procedure may allow a density reconstruction which is completely 
independent of any other method that has thus far been proposed. We have demonstrated its 
effectiveness in simulated spin catalogs. Due to the nature of this letter, only the key results 
essential to the reconstruction algorithm have been presented. Detailed intermediate steps and 
justifications will be shown elsewhere (Lee k Pen 2000). 

We thank Giuseppe Tormen and Antonaldo Diaferio for providing the halo-finding algorithm, 
and the referee for useful suggestions. J. Lee also thanks Anatoly Klypin for helpful comments 
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